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(57) Abstract 

A method approximates propagation delay through a logic device. Operation of the logic device is divided into a first re- 
gion and a second region. A boundary between the first region and the second is based on duration of input ramp to the logic de- 
vice and amount of capacitive load driven by the logic device. For example, the boundary between the first region and the second 
occurs where for each value of the capacitive load, an output ramp for the logic device is one hair complete when the input ramp 
is complete. When the logic device operates in the first region, a first formula is used to obtain a first value representing delay 
through the logic device. The first formula varies the first value based on the duration of the input ramp to the logic device and 
the capacitive load driven by the logic device. When the logic device operates in the second region, a second formula is used to 
obtain the first value. The second formula also varies the first value based on the duration of the input ramp to the logic device 
and the capacitive load driven by ihe logic device. 
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TIMING MODEL AND CHARACTERIZATION SYSTEM FOR LOGIC 
SIMULATION OF INTEGRATED CIRCUITS 

Technical Field 

The present invention concerns an accurate timing model for logic 
5 simulation of integrated circuits. The timing model reflects the timing 
delays caused by input ramp and output load of circuitry. 

Background Art 

When designing an integrated circuit, it is generally desirable to 
simulate the functioning of logic circuitry within the integrated circuit. In 

10 order to accurately access the performance of the circuitry, it is additionally 
desirable that the logic simulation of logic include an accurate assessment 
of timing delays through the circuitry, at least in the critical paths. 

In general, timing delays through circuitry are caused by propagation 
delays through and between logic cells which comprise the circuitry. The 

15 actual amount of propagation delay through and between logic cells is 

generally dependent on various capacitances within and between the logic 
cells, as well as the current available to charge or discharge the 
capacitances. 

When developing a timing model, it is desirable to take into account 
20 the time delay introduced by charging the input capacitance of logic cells. 
This time delay is directly affected by the input current available to charge 
the input capacitance to the logic cell. This input current, in turn, is 
directly affected by the fan out of the output of the logic cell providing the 
input current. 

25 In the prior art, timing delay caused by charging input capacitance 

was approximated as being directly proportional to the output fan out of the 
logic cell providing current for the input. However, this has proved to be an 
. inaccurate modelling technique because propagation delay through the logic 



cell is additionally influenced by the shape of the voltage input ramp of the 
input signal to the logic cell. The effect of the shape of the voltage input 
ramp on the timing delay may be approximated by using general piecewise 
linear functions of the voltage input ramp and the output load. However, 
5 this approach does not easily lead to a practical characterization system. In 
order to obtain a fairly accurate timing model, a circuit simulation needs to 
be performed for every grid point of the function. A sufficient number of grid 
points need to selected to obtain the desired accuracy of the timing model. 
This can require numerous circuit simulations. 

10 

In accordance with the preferred embodiment of the present ' 
invention, a method is presented for approximating propagation delay 
through a logic device. Operation of the logic device is divided into a first 

15 region and a second region. A boundary between the first region and the 
second is based on duration of input ramp to the logic device and amount of 
capacitive load driven by the logic device. For example, the boundary 
between the first region and the second occurs where for each value of the 
capacitive load, an output ramp for the logic device is one half complete 

20 when the input ramp is complete. 

When the logic device operates in the first region, a first formula is 
used to obtain a first value representing delay through the logic device. The 
first formula varies the first value based on the duration of the input ramp to 
the logic device and the capacitive load driven by the logic device. When the 

25 logic device operates in the second region, a second formula is used to obtain 
the first value. The second formula also varies the first value based on the 
duration of the input ramp to the logic device and the capacitive load driven 
by the logic device. 



3 

For example, the first formula may have a form as follows: 
DTo = Aq + B*(INPUT RAMP) + Do*(LOAD) 
where DTo represents the first value, INPUT RAMP represents duration of 
the input ramp, LOAD represents capacitive load driven by the logic device, 
5 and An, B and Do are coefficients determined using circuit simulations. 
Similarly, the second formula may have a form as follows: 
DTo = Ao+B*(CRITINRAMP)+Do*(LOAD)+Z*(INPUT RAMP - 

CRITINRAMP) 

where CRITINRAMP represents INPUT RAMP at the border between the 
10 first region and the second region for the LOAD, and Z is determined using 
circuit simulations. 

The above examples for the first formula and the second formula use 
first order linear approximations for the input ramp of the logic device. 
Alternately a second order approximation of the input ramp may include 
15 "tails" of the input ramp. Using the second order model, the first formula 
' may have a form as follows: 

DTo = Ao + FAT*B*(INPUT RAMP) + Do*(LOAD) 
wherein DTo represents the first value, INPUT RAMP represents duration 
of the input ramp, LOAD represents capacitive load driven by the logic 
20 device, and Ao, B, Do and FAT are coefficients determined using circuit 
simulations. 

Similarly, using the second order model, the second order model may 

'*■-,- 

have a form as follows: 

DTo = DTo(CRITINRAMP) '+ Z * (INPUT RAMP - CRITINRAMP) 
25 wherein CRITINRAMP represents INPUT RAMP at the border between the 
first region and the second region for the LOAD, Z is a coefficient 
determined using circuit simulations, and DT(CRITINRAMP) is the first 



value when the logic ramp has an input ramp equal to CRITINRAMP and a 
load equal to LOAD. 

In the preferred embodiment of the present invention, the method also 
includes calculating a next input ramp for a next logic device coupled to an 
5 output of the logic device. When the logic device operates in the first region, 
a third formula is used to generate a second value representing duration 
from beginning of the input ramp of the logic device to an end of an output 
ramp of the logic device. The second value is based on the length of the input 
lamp to the logic device and the amount of the load driven by the logic device. 
10 When^the logic device operates in the second region, a fourth formula is 
used to generate the second value. The first value is then subtracted from 
^ the serond value to ob • " ? ^ - 

For example, the third formula may have a form as follows: 
DTi = Ai + B*(INFUT EAMP) + D]*(LOAD) 
15 where DTi represents the second value, and Ai and Di are coefficients 
determined using circuit simulations. Similarly, the fourth formula may 
have a form as follows: 

DTi = Ai + B*(INPUT RAMP) + DrKLOAD). 
Alternately, using the second order model which takes into account 
20 "tails" of the input ramp, the third formula may have a form as follows: 
DTi = Ai + FAT*B*(INPUT RAMP) + Di*(LOAD) 
wherein DTi represents the second value, and Ai and Di are coefficients 
determined using circuit simulations. Similarly, using the second order 
model, fourth formula may have a form as follows: 
25 DTi = DTi(CRITINRAMP) + B * (INPUT RAMP - CRITINRAMP). 

The present invention is advantageous over the prior art in that it 
allows for a fairly accurate timing model of a logic device with the need to 
perform only a limited number-of circuit simulations. 
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Brief Pescriptign of flift Prawings 

Figure 1 shows an example of circuitry through which propagation 
delay may be calculated. 
5 Figure 2 shows a graph on which is superimposed a ramp waveform. 

Figure 3 shows a graph on which is superimposed a ramp waveform. 
Figure 4 shows two operation regions for a logic device in accordance 
with the preferred embodiment of the present invention. 

Figure 5 shows a graph on which is superimposed an input ramp 
10 waveform and a resulting output ramp waveform in a first operation region 
shown in Figure 4 in accordance with the preferred embodiment of the 
... present invention. .,...^,..4, ,, r .- v ,v ^<^ r .*v^*-; -^.«-p>y*?*. ^^mv^y^®-* ■ ^•-••42>v--.-. '^-^ 
Figure 6 shows a graph on which is superimposed an input ramp 
waveform and a resulting output ramp waveform in a second operation 
15 region shown in Figure 4 in accordance with the preferred embodiment of 
the present invention. 

Figure 7 shows a circuit used for a second order model simulation in 
accordance with the preferred embodiment of the present invention. 

20 Description of the Preferred EmMangnt 

There are many causes of propagation delay through logic circuitry. 
For example, Figure 1 shows circuitry which includes a logic device 11, a 
logic device 12, a logic device 13, a logic device 14, a logic device 15, a logic 
device 16, a logic device 17, a logic device 18, a logic device -19, a logic device 
25 20 and a logic device 21 connected as shown. For simplicity in the discussion 
of the present invention, inverters as used for the logic devices. However, 
the present invention is applicable to all logic devices used in the 
manufacture of integrated circuits. 



Causes for propagation delay through the circuitry shown in Figure 1 
include, for example, the time required to charge or discharge internal 
capacitance within each of logic devices 11-21, the time required to charge or 
discharge input capacitance to each of logic devices 11-21, and the time 

5 required to overcome inductance in the lines hetween logic devices 11-21. 
The present invention, however, focuses on the effects that input transition 
time and output load has on propagation delay. The input transition time of 
a logic device is the time required to charge or discharge input capacitance 
of the logic device. The output load is the total amount of capacitance which 

10 the logic device has to charge or discharge. The output load is directly 
affected by, for example, the fan out of the output of the logic device. 
- ••••-v For example, as shown in Figure 1; logic device 11 -has a fan out of - 
two. Logic device 12 has a fan out of three. Logic devices 13-17 have a fan out 
of one- Because of the relatively large fan out of logic device 12, it is expected 

15 that the output transition time of (and thus the propagation delay through) 
logic device 12 will take longer than, for example, the output transition time 
of logic device 13, which has a fan-out of only one. 

Figure 2 shows a graph 129 on which is superimposed a ramp 
waveform 153 which represents a voltage ramp during a transition period. 

20 A vertical axis 128 of graph 129 represents voltage at an input or an output of 
a logic device. A horizontal axis 127 represents passage of time. A zero 
point 130 of graph 129 indicates zero voltage and the passage of zero 
nanoseconds. A point 141 on horizontal axis 127 represents the passage of 
2.5 nanoseconds. A point 142 on horizontal axis 127 represents the passage 

25 of 5.0 nanoseconds. A point 143 on horizontal axis 127 represents the 

passage of 7.5 nanoseconds. A point 144 on horizontal axis 127 represents 
the passage of 10.0 nanoseconds. A point 131 on vertical axis 128 represents 
0.5 volts. A point 132 on vertical axis 128 represents 1.0 volts. A point 133 on 
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. vertical axis 128 represents 1.5 volts. A point 134 on vertical axis 128 
represents 2.0 volts. A point 135 on vertical axis 128 represents 2.5 volts, A 
point 136 on vertical axis 128 represents 3.0 volts. A point 137 on vertical axis 
128 represents 3.5 volts. A point 138 on vertical axis 128 represents 4,0 volts. 

5 A point 139 on vertical axis 128 represents 4.5 volts. A point 140 on vertical 
axis 128 represents 5.0 volts. 

Ramp waveform 153 shows voltage verses time for a transition from 
low to high. Ramp waveform 153 may be approximated by a waveform 154, 
Waveform 154 approximates as a straight line the transition of ramp 

10 waveform 153 from a low voltage to a high voltage. The total transition time 
155, or ramp, of ramp waveform 153 is defined as the duration of time 
between a time 151 (To), which is the point at which waveform 154 begins to 
increase in voltage, and a time 152 (Tx), which is the point at which 
waveform 154 ceases to increase in voltage. 

15 Figure 3 shows a ramp waveform 163 superimposed on graph 129. 

Ramp waveform 163 shows voltage verses time for a transition from high to 
low. Ramp waveform 163 may be approximated by a waveform 164. 
Waveform 164 approximates as a straight line the transition of ramp 
waveform 163 from a high voltage to a low voltage. The total transition time 

20 165, or ramp, of ramp waveform 163 is defined as the duration of time 

between a time 161 (To), which is the point at which waveform 164 begins to 
decrease in voltage, and a time 162 (Ti), which is the point at which 
waveform 164 ceases to decrease in voltage. 

The terminology introduced in the discussion of Figures 2 and 3 may 

25 be used in the general description of input and output transition timing for a 
logic cell. For an input pin i (not shown) of a logic cell, the beginning of a 
ramp is defined as To,i and the end of the ramp is defined as Ti,i. For an 
output pin j (not shown) of a logic cell, the beginning of a ramp is defined as 



To j and the end of the ramp is defined as Ti j. Using this terminology, 50%- 
50% propagation delay (Tau) between input pin i and output pin j, is given by 
the following Formula 1: 

Formula 1 

5 Tau,ij = l/2(Ti j + To j) - + T 0,0 

While it is often the case that propagation delay is calculated using 

Tau, for the present invention, other calculated values may also be used. 

For example, the value DTo is defined as the difference in time between the 

beginning of the ramp for the input and the beginning of the ramp for the 
10 output^ More generally, the value DTo between input pin i and output pin j 

is described by the following Formula 2: 

' DTo,ij=T 0 j-To,i 
The value DTi is defined as the difference in time between the beginning of 
15 the ramp for the input and the end of the ramp for the output. More 

generally, the value DTi between input pin i and output pin j is described by 
the following Formula 3: 

Formula 3 
DTi4j=Tij-T 0 ,i 

20 For the present invention, Wo regions of operation of a logic device are 

considered. These regions are shown in Figure 4. In Figure 4, a vertical 
axis 171 represents increases in output load for the logic device. An axis 172 
represents increase in time duration for the input ramp of the input ramp. 
A first region 174 is referred to as the fast input ramp region because in this 

25 region the input ramp rises relatively fast in comparison to the output ramp. 
A second region 175 is referred to as the slow input ramp region because in 
this region the input ramp rises relatively slowly in comparison to the 
output ramp, A dividing line 173 between fast input ramp region 174 and 
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slow input ramp region 175 may be determined empirically as further 
discussed below. 

Figure 5 shows an example of an input ramp 251 and a resulting 
output signal waveform 252 in fast input ramp regions 174. As may be seen 

5 from Figure 5, input ramp 251 falls at a relatively fast rate as compared to 
the rate of the rise of output signal waveform 252. 

Input ramp 251 and output signal waveform 252 are superimposed on 
a graph 229. A vertical axis 228 of graph 229 represents voltage. A 
horizontal axis 227 represents passage of time. A zero point 230 of graph 229 

10 indicates zero voltage and the passage of zero nanoseconds. A point 241 on 
horizontal axis 227 represents the passage of 1.67 nanoseconds. A point 242 
on horizontal axis 227 represents the passage of 3.33 nanoseconds. A point 
243 on horizontal axis 227 represents the passage of 5.00 nanoseconds. A 
point 244 on horizontal axis 227 represents the passage of 6.67 nanoseconds. 

15 A point 245 on horizontal axis 227 represents the passage of 8.33 

nanoseconds. A point 246 on horizontal axis 227 represents the passage of 
10.0 nanoseconds. A poini 231 on vertical axis 228 represents 0.5 volts. A 
point 232 on vertical axis 228 represents 1.0 volts. A point 233 oh vertical axis 
228 represents 1.5 volts. A point 234 on vertical axis 228 represents 2.0 volts. 

20 A point 235 on vertical axis 228 represents 2.5 volts. A point 236 on vertical 
axis 228 represents 3.0 volts. A point 237 on vertical axis 228 represents 3.5 
volts. A point 238 on vertical axis 228 represents 4.0 volts. A point 239 on 
vertical axis 228 represents 4.5 volts. A point 240 on vertical axis 228 
represents 5.0 volts. 

25 Figure 6 shows an example of an input ramp 281 and a resulting 

output signal waveform 282 in slow input ramp regions 175. As may be seen 
from Figure 6, input ramp 281 falls at a relatively slow rate as compared to 
the rate of the rise of signal waveform 282. 



Input ramp 28 1 and output signal waveform 282 are superimposed on 
a graph. 259. A vertical axis 258 of graph 259 represents voltage. A 
horizontal axis 257 represents passage of time. A zero point 260 of graph 259 
indicates zero voltage and the passage of zero nanoseconds. A point 271 on 

5 horizontal axis 257 represents the passage of 1.25 nanoseconds. A point 272 
on horizontal axis 257 represents the passage of 2,5 nanoseconds. A point 
273 on horizontal axis 257 represents the passage of 3.75 nanoseconds. A 
point 274 on horizontal axis 257 represents the passage of 5.0 nanoseconds. 
A point 275 on horizontal axis 257 represents the passage of 6.25 

10 nanoseconds. A point 276 on horizontal axis 257 represents the passage of 
7.5 nanoseconds. A point 261 on vertical axis 258 represents 0,5 volts. A 
point 262 on vertical axis 258 represents volts. A prnnt 263 on vertical axis 
258 represents 1.5 volts. A point 264 on vertical axis 258 represents 2.0 volts. 
A point 265 on vertical axis 258 represents 2.5 volts. A point 266 on vertical 

15 axis 258 represents 3.0 volts. A point 267 on vertical axis 258 represents 3.5 
volts. A point 268 on vertical axis 258 represents 4.0 volts. A point 269 on 
vertical axis 258 represents 4.5 volts. A point 270 on vertical axis 258 
represents 5.0 volts. 

In fast input ramp region 174, input ramp 251 is completed switched 

20 before output signal waveform 252 significantly changes. For a CMOS 

technology implementation of an inverter, one transistor of the output buffer 
of the inverter would be on in saturation and the other transistor of the 
output buffer of the inverter would be off during substantially all the time 
output signal waveform 252 switches. This is essentially equivalent to a 

25 constant current source charging or discharging the load capacitance. In 
the preferred embodiment of the present invention, for fast input region 174, 
value DTo and value DTi are approximated using the following Formula 4 
and Formula 5, respectively: * 
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Fflrmylfl 4 

DTo = Ao + B*(INPUT RAMP) + Do*(LOAD) 
Fprmydfi 5 

5 DTi = Ai + B*(INPUT RAMP) + Di*(LOAD) 

In Formulas 4 and 5, values Ao, Ai, B, Do and Di are coefficients which 
may be determined empirically. 

In slow input ramp region 175, input ramp 251 switches so slowly that 
output signal waveform 252 switches while input ramp 252 is still switching. 

10 For a CMOS technology implementation of an inverter, both transistors of 
the output buffer of the inverter are on when output signal waveform 252 
switches. In the preferred embodiment of the present invention, for slow 
input region 175, value DTo and value DTi are approximated using the 
following Formula 6 and Formula 7, respectively: 

15 Fprm\4a 6 

DTo = Ao+B*(CRITINRAMP)+Do*(L6ADHZ*(INPUT RAMP - 

CRITINRAMP) . 

Fonrndfl 7 

20 DTi = Ai + B*(INPUTRAMP) + Di*aOAD) 

In formula 6, the additional co-efficient Z is a proportionality factor 

, calculated from the DC transfer function between the input voltage and the 
output voltage. In the slow ramp region, a plot of DTo versus the input ramp 
is a straight line for a given output load. Z is the slope of that line and can be 

25 calculated from the DC transfer curve. CRITINRAMP is the INPUT RAMP 
for each load which is along dividing line 173 between fast input ramp 
region 174 and slow input ramp region 175. In the preferred embodiment of 
the present invention, the CRITINRAMP (and thus dividing line 173) is 
defined to be at the point for each load where the end of the input ramp 

30 occurs at 50% of the output ramp. Since time 0 is considered to be at the 
beginning of the input ramp, then the end of the input ramp occurs at 



. INPUT RAMP. Therefore, CRITINRAMP occurs when INPUT RAMP 
equals 1/2 * (DTo + DTlX CRITINRAMP is therefore defined by the 
following Formula 8: 

Formula 8 

5 CRrriNRAMP = l/2*(DTo + DTi) 

The coefficients for the formulas may be calculated using circuit 
simulations. First, three circuit simulations are performed to determine 
values for the coefficients (Ao, Ai, B, Do and Di) used in formula 4 and 
formula 5. In the circuit simulations values for load and input ramp are 
10 chosenso that the three circuit simulations will all fell within fast input 
ramp region 174. Actually, initially the location of dividing line 173 is not 
known so that the actual boundaries of fast input ramp region 174 are not 
known. Therefore, later in the model it is necessary to perform a check to 
see if the initial three circuit simulations were in fast input ramp region 
15 174. If not, the circuit simulations are performed again using a faster input 
. ramps and/or higher loads. 

For each simulation, using formulas 2 and 3, DTo and DTi are found. 
Then using formulas 4 and 5, Ao, Ai, B, Do and Di may be calculated A 
first consistency check may be done by solving formulas 4 and 5 separately 
20 for the coefficient B. If the coefficient B does not have the same value, then at 
least one of the circuit simulations was not in fast input ramp region 174. 

From the calculated values for coefficients Ao, Al, B, Do and Di, 
using Formula 8, CRITINRAMP is calculated. Once CRITINRAMP is 
found, the chosen circuit simulations may again be checked to be sure they 
25 fall within fast input ramp region 174. 

An additional circuit simulation is performed to find additional points 
on the DC transfer curve. From these points, the DC transfer function may 
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be approximated. A linear approximation of the DC transfer function yields 
the proportionality constant Z. 

In a manner similar to that described above, coefficients are 
calculated for each logic device type in a logic device library. Delay through 

5 a circuit built from logic devices in the logic device library may then be 
calculated. Total propagation delay through a circuit is the total amount of 
time for a signal to propagate through the circuit. The path of logic cells in a 
circuit through which a signal takes longest to propagate is called the 
critical path of the circuit. Total propagation for a circuit may be calculated 

10 by summing up the total delays through all logic cells in the critical path of 
the circuit. 

Formula 1, for example, may be used to calculate delay through each 
logic cell. Additionally 9 starting with the input of the circuit and continuing 
through to the output of the circuit, the input ramp and the resulting output 

15 ramp (DTi - DTo) may be calculated for each logic device in the logic path. 
The output ramp for one logic device is the input ramp for the next logic 
device. The output ramps are separately calculated for both rising inputs 
and falling inputs. 

To more fully illustrate how the present invention is practiced, an 

20 example is given below for the case of an inverter. In the example, three 
simulations are done with the input ramp and output load (LOAD) listed in 
Table 1 below: 



INPUT HAMP LO AD 

25 .25 nS lpF 

.25 nS 2pP 

InS 2pP 



For the first two simulations input voltage (VTNPUT) is given by the 
following Formula 9: 



Formula 9 

VINPUT = 18.2 * TIME For TIME < .25 nS 

VTNPUT = 4.55 For TIME > .25 nS 

Similarly, for the third simulation, input voltage (VINPUT) is given by the 

5 following Formula 10: 

Formula 10 

VINPUT = 4.55* TIME For TIME < 1 nS 

VINPUT = 4.55 For TIME >lnS 

In the above Formula 9 and Formula 10, TIME is measured from the 
10 beginning of the INPUT RAMP and VINPUT is measured in volts. 

For each of the simulations, a linear approximation is made of the 
output of the simulated circuit. The linearization is done by a two point fit 
through times at which the output voltage is at 3.15 volts and at 1.4 volts. In 
the simulation, the results of the linearization are given by Table 2 below: 

15 Table 2 

Time when Time when 

INPUT RAMP LOAD OUTPUT at 3.15 V OUTPUT = 1.4 V 

.25 nS IpF L027nS 2.107nS 

-25 nS 2pF L806nS 3.911nS 

20 InS 2pF 2.289 nS 4.394nS 

Using linear approximation for the first simulation gives an output 

voltage (VOUTPUT) for the first simulation as given by Formula 11 below: 

Formula 11 

VOUTPUT = 4.55 For TIME < .163 nS 

25 VOUTPUT = 4.55 - 4.55 * (TIME - .163)/ (2.971 - .163) 

For .163 nS < TIME < 2.971 nS 
VOUTPUT = 0 For TIME > 2.971 nS 

Using linear approximation for the second simulation gives an output 

voltage (VOUTPUT) for the second simulation as given by Formula 12 below: 

30 

Formula 12 

VOUTPUT = 4.55 For TIME < .122 nS 

VOUTPUT = 4.55 - 4.55 * (TIME - .122)/ (5.595 - .122) 

For .122 nS < TIME < 5.595 nS 
35 VOUTPUT = 0 For TIME > 5.595 nS 



WO 93/18468 ,'€-Y" - PCT/US93/01747 

15 

Using linear approximation for the third simulation gives an output voltage 
(VOUTPUT) for the third simulation as given by Formula 13 below: 

Formula 13 

VOUTPUT = 4.55 For TIME < 1.251 nS 

5 VOUTPUT = 4.55 -4.55* (TIME -1.251)/ (6.724 -1.251) 

For 1.251 nS < TIME < 6.724 nS 
VOUTPUT = 0 For TIME > 6.724 nS 

Values for DTo and DTi for the three simulations are extracted from the 
linear approximations as listed in Table 3 below: 



10 Table 3 

INPUT RAMP LOAD DT£ XHl 

.25 nS lpF .163 nS 2.971nS 

.2frnS 2pF .122 nS 5.595 nS 

InS 2pF 1.251nS 6.724 nS 



15 Now the values for DTo and DTi may be substituted into Formula 4 

and Formula 5 to solve for coefficients An, Ai, B, Do and Di. For example, 

substituting the values for INPUT RAMP, LOAD and DTo for the three 

simulations into Formula 4 gives the following equations: 

.163 nS = Ao + B*U5 nS) + Do*(l pF) 
20 .122nS = Ao+B*(.25nS) + Do*(2pF) 

1.251 nS = Ao + B*(1,0 nS) + Do*(2 pF) 

These three equations are solved for Ao, B, and Do to produce the following 

values: 

Ao=043 

25 B = .644 

Do = -.041 

Similarly, substituting the values for INPUT RAMP, LOAD and DTi for the 
three simulations into Formula 5 gives the following equations: 

2.971 nS = Ai + B*(.25 nS) + Dl*(l pF) 
30 5.595 nS = Ai + B*(.25nS) + Di*(2pF) 

6.724 nS = Ai + B*(1.0nS) + Di*(2 pF) 



These three equations are solved for Ai, B, and Di to produce the following 
values: 

Ai = .186 
B = .644 

5 Di = 2.624 

The first consistency check is done by noting that the coefficient B has the 
same value for the solutions of both Formula 4 and Formula 5. If two 
different values of B had been obtained, it would imply that at least one of the 
original simulations was not in the fast ramp region. The simulations 

10 would have to be redone with faster input ramps and/or larger output loads. 

'Now that the coefficients have been determined CRITINRAMP may be 

calculated from Formula 8. For example, Formula 8 is simplified as below 

(remember that at CRTFINRAMP: INPUT RAMP = CRITINRAMP): 

CRTFINBAMP = 1/2 * (DTo + DTi) 
15 CRITINRAMP = 1/2* DTo + 1/2 *DTi 

CRITINRAMP = 1/2 * [Ao + B*(INPUT RAMP) + Dn*(LOAD)] + 1/2 * [Ai + 

B*(INPUT RAMP) + Di*(LOAD)] 
CRITINRAMP = 1/2 * (Ao + Ai) + 1/2 * (Do + Di)*(LOAD) + 
B*(CRITINRAMP) 
20 (mmNRAMP = l/2*[Ao + Ai + (Do + Di)*(LOAD)]/(l-B) 

Using the values for the coefficients Ao, Ai, B, Do and Di determined 

above, CRnTNRAMP has the value given in Formula 14 below: 

Formula 14 
CRITINRAMP = .322 + 3.628 * (LOAD) 

25 The second consistency check is done by checking to see Whether the three 

simulations fall within fast input ramp region 174. For the two values for 

loads used in the three simulations, CRITINRAMP is given by Table 4 

below: 

Table 4 

30 LOAD CRITINRAMP 

lpF 3.95QnS 
2pF 7.578 nS 
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The first simulation was done with a 1 pF load and an input ramp of .25 nS. 
Table 4 indicates this is within fast input ramp region 174. The second and 
third simulations were done with a 2 pF load and input ramp of .25 nS and 1 
nS, respectively. Table 4 indicates these are also within fast input ramp 
5 region 174. 

An additional circuit simulation of DC transfer curves and linear 
approximations is made to obtain the proportionality constant Z. In the 
circuit simulation, DC input voltages are found which cause output voltage 
to be 3.15 volts and 1.4 volts. The circuit simulation yields the values given in 
10 Table-5 below: 

Table 5 

Output, Voltage ftnwt Voltage 

3.15 V 2.176 V 

1.40 V 2.298 V 

15 Using linear approximations, the DC transfer curve is calculated as 

given in Formula 15 below: 

Forpiyte 15 

VOUTPUT = 4.55 For VINPUT < 1.251 

VOUTPUT = 4.55 - 4.55 * (VINPUT - 2.078Y (2.396 - 2.078) 
20 For 2.078 < VINPUT < 2.396 

VOUTPUT = 0 For TIME > 6.724 nS 

In the slow ramp region, a plot of DTo versus the input ramp is a 

straight line for a given output load. Z is the slope of that line and can be 

calculated from the DC transfer curve. In determining proportionality 

25 constant Z it is helpful to define the following breakpoints. Breakpoint V + is 
the value for VINPUT at which VOUTPUT reaches 4.55. Breakpoint V" is 
the value for VINPUT at which VOUTPUT reaches 0. For the present 
example, from Formula 15, it is seen that breakpoint V+ is 2.078 volts and 
breakpoint V' is 2.396. Thus, proportionality constant Z may be calculated 

30 from Table 6 below: 



Table 6 

VINPUT VOUTPUT Z 

RISING RISING (V--0)/ (4.55-0) 

RISING FALLING (V+ - 0) / (4.55 - 0) 

5 FALLING RISING (4.55 -V")/ (4.55-0) 

FALLING FALLING (4.55 - V+) / (4.55 - 0) 

In the given example, the logic device is an inverter with input rising. 

Therefore, VINPUT is rising and VOUTPUT is falling. From Table 6, the 

value of proportionality constant is therefore; 

10 : 2 = 2.078/4.55 = ,457 

In the logic device library for the logic device, the values for coefficients Aq, 
Al, Brt>0 and Di are stored. Additionally, the value for CRfflNRAMP from 
Formula 14 is stored and four values calculate from Table 6 are stored for 
proportionality constat Z. When calculating propagation delay through a 

15 first logic device, DTo of the first logic device is equal to the beginning of an 

input ramp for a second logic device following the first logic device. 
Similarly the value (DT^ - DTo) f° r the logic device is equal to the input 
ramp for the second logic device. By determining DTo ^ DTi for each logic 
device and by using the value (DT^ - DTo) as 311 input ramp for a next logic 

20 device, the total propagation delay through circuitry is accomplished 

In the above example, the linear approximations assumed perfectly 
linear inputs. However, tails in the input waveform were not accounted for. 
The tail in the input waveform has a noticeable effect in fast input region 
174; however, the tail in the input waveform has a negligible effect in slow 

25 input ramp region 174. To take account of this second order effect, two 
factors to adjust for the tails (FAT) may be included, one for rising inputs 
and one for falling inputs. To take these into account, Formula 4 and 
Formula 5 are modified resulting, respectively, in Formula 16 and Formula 
17 given below: 
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i. 

Formula Ifi 

DTq = Ao + FAT*B*(INFUT RAMP) + Do*(LOAD) 
Formula 17 

5 DTi=Ai + FAT*B*(INPUTRAMP) + Di*(LOAD) 

As in the above example, the value FAT may be determined by using circuit 
simulations. Figure 7 shows a model used for circuit simulations for the 
second order model. Timing delay is modeled for an inverter 30. A 
capacitor 36 connected between an output 35 of inverter 30 and a ground 37 

10 represents output load of inverter 30. A capacitor 34 connected between an 
input 33 of inverter 30 and ground 37 represents input load of inverter 30, 
An inverter 32 is used in order to place on input 33 of inverter 30 a realistic 
voltage, that is an input voltage with a tail. 

Three simulations are performed using a falling (linear) input on an 

15 input 31 of inverter 32 with a ramp of .1 nS. The simulations are performed 
as the circuit simulations. The simulations are performed as described in 
the prior examples, with the exception each of the three simulations 
includes a simulation of both inverter 32 and inverter 30. In a first 
simulation, capacitance 34 (Cint) has a value of 0 pF and capacitance 35 

20 (Cout) has a value of 1 pF. The first simulation of inverter 32 and linear 
approximation yields the following voltage (Vint) at input 33 of inverter 30, 
shown by Formula 18 below: 

Fgrmylfl 18 

Vint = 0 For TIME <. 064 nS 

25 Vint = 4.55* (TIME -.064)7 (.395 -.064) 

For .064 nS < TIME < .395 nS 
Vint = 4.55 For TIME >. 395 nS 

The first simulation of inverter 30 and linear approximation yields the 

following voltage (Vout) at output 35 of inverter 30, shown by Formula 19 

30 below: 



Formula 19 

Vout = 4.55 For TIME < ,311 nS 

Vout = 4.55 - 4.55 * (TIME - .311)/ (3.126 - .311) 

For .311 nS < TIME < 3.126 nS 
5 Vout = 0 For TIME > 3.126 nS 

In Formula 18 and Formula 19, TIME is measured from the beginning of 

the ramp on input 31 of inverter 32. Inverter 30, therefore, on input 33 sees 

an input transition that starts at .064 nS and ends at .395 nS. The output 

transition on output 35 of inverter 30 starts at .311 nS and ends at 3.126 nS. 

10 Thus the input ramp on input 33 of inverter 30 is .395 minus .064, that is .331 

nS. Using the Formula 2, DTo is seen to be .311 minus .064, that is, .247 nS. 

Using Formula 3, DTi is seen to be 3.126 minus .064, that is, 3.062 nS. 

Table 7 below summarizes the results of all three simulations: 

15 £mt £fi3i£ INPUT RAMP DJfr Ull 

OpF lpF .331nS .247 nS 3.062 nS 

OpF 2pF .331 nS .206 nS 5.686 nS 

.2pF 2pF .826 nS .826 nS 6.040 nS 

20 The above values are used to solve for the coefficients ao, ai, b, do and 

di in the following following Formula 20 and Formula 21: 

Formula 20 
DTo = ao + b*ONPUT RAMP) + do*(Cout) 

25 Formula 21 

DTi = di + b*(INPUT RAMP) + di*(Cout) 

Formula 20 and Formula 21 uses the coefficients ao, ai, b, do and di because 
the values of the coefficients sometimes vary from the values of the 
30 coefficients Ao, Ai, B, Do and Di used in the first order model. As above, the 
values in Table 7 provide for three sets of equations for Formula 20 and 
Formula 21. Solving the three equations yields the following values: 
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ao = .051 
ai = .201 
b = .715 
do = -.041 

5 di = 2.624 

When these' are compared with the values for Ao, Ai, B, Do and Di 
calculated above, it is seen that the values for do and di are equivalent to the 
values for Do and Di, respectively. The values for ao and ai vary slightly 
from the values for Ao and Ai, respectively. The most significant difference 

10 is between the value for b and the value for B. The value FAT is obtained by 
the ratio of b to B, as shown by Formula 22 below: 

Fgrmylfl 32 

FAT =b/B 

= .715/. 644 

15 -•■< - / =111^^^^^ " . >: ' — ; " ' 

Once calculated, the value for FAT is stored in the logic device library along 
with the values for the coefficients Ao, Ai, B, Do and Di. The value for FAT 
is dependent on the process technology and does not depend on the logic 
design of a logic cell. Therefore, a single entry for the value of FAT may be 

20 obtained from characterization of a single logic cell and stored for use for all 
logic cells in a cell library. 

The value for CRITINEAMP for the second order model may be 
recalculated from Formula 8, to be as set out in Formula 23 below: 

FpiTnylfl 23 

25 CRITINEAMP = 1/2* [Ao + Ai + (Do .+ DiMLOAD)] / (1 - FAT*B) 

Using the values for the coefficients Ao, Ai, B, Do, Dl and FAT 
determined above, CRITINRAMP has the value given in Formula 24 below: 

Egnbala 24 

CRITINRAMP = .402 + 4.532 * (LOAD) 
30 Using the second order model, the following Formula 25 and Formula 26 are 
used to calculate DTo and DTi in slow input ramp region 175: 



Formula 25 

DTo = DTo(CRITINRAMP) + Z * (INPUT RAMP - CRITINRAMP) 

Formula 26 

5 DTi .= DTi(CRITINRAMP) + B * (INPUT RAMP - CRITINRAMP) 

Li Formula 25, DTo(CRITINRAMP) is the value of DTo at CRITINRAMP. 

In Formula 26, DTi(CRITINRAMP) is the value of DTi at CRITINRAMP. 
Once the logic device has been characterized, as above, it is possible to 

determine a value for DTo and DTi for any input ramp and load. For 
10 example, suppose the load (LOAD) is equal to .5 pF and the input ramp is 

equal to .25 nS. Then, using Formula 24, CRITINRAMP has the following 

r- 

value: 

CRITINRAMP = .402 + 4.532 * (LOAD) 
CRITINRAMP = .402 + 4.532 * (.5) 
15 CRITINRAMP = 2.668 nS 

Since the input ramp (.25 nS) is smaller than the CRITINRAMP 

(2.668 nS), the example falls in fast input ramp region 174. Therefore, 

Formula 16 is used to calculate DTo and Formula 17 is used to calculate 

DTi, as shown below: 

20 DTo = Ao + FAT*B*(INPUT RAMP) + Do*(LOAD) 

DTo = .043 + 1.11 * .644 * (-25) - -041 * (.5) 
DTo = .201nS 

DTi = Ai + FAT*B*(INPUT RAMP) + Di*(LOAD) 
25 DTi = .186 + 1.11 * .644 * (.25) + 2.624 * (.5) 

DTi = 1.677 nS 

As an alternate example, suppose the load (LOAD) is equal to .5 pF and the 
input ramp is equal to 3 nS. Then CRITINRAMP still has the following 
value: 

30 CRITINRAMP = .402 + 4.532 * (LOAD) 

CRITINRAMP = .402 + 4.532 * (.5) 
CRITINRAMP = 2.668 nS 

Since the input ramp (3 nS) is now larger than the CRITINRAMP 

(2.668 nS), the example falls in slow input ramp region 175. Therefore, 
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Formula 25 is used to calculate DTo and Formula 26 is used to calculate 
DTi. First DTo (CRITINRAMP) and DTi(CRITINRAMP) are calculated as 
shown below: 

DTo(CRITINRAMP) = Ao + FAT*B*(CRITINRAMP) + Do*(LOAD) 
5 DTo(CRITINRAMP) = 043 + 1.11*.644*(2.668) - .041 * (.5) 

DTo(CRITINRAMP) = 1930 nS 

DTi(CRITINRAMP) = Ai + FAT*B*(CRITINRAMP) + Di*(LOAD) 
DTi(CRITINRAMP) = -186 + 1.11*.644*(2.668) + 2.624 * (.5) 
10 DTi(CRITINRAMP) = 3.406 nS 

Once PTo(CRITINRAMP) and DTi(CRITINRAMP) have been calculated, 

Formula 25 is used to calculate DTo and Formula 26 is used to calculate 

DTl, as shown below: 

DTo = DTo(CRITINRAMP) + Z * (INPUT RAMP - CRITINRAMP) 
15 > DTo = 1-930 + .457 * (3 - 2.668) : - • . • 

DTo = 2.082 nS 

DTi = DTKCRITINRAMP) + B * (INPUT RAMP - CRTTINRAMP) 
DTl = 3.406 + .644 * (3 - 2.668) 
20 DTi = 3.620 nS 

As discussed above the input ramp for the following logic device is equal to 
DTi -DTo- 

The foregoing discussion discloses and describes merely exemplary 
' methods and embodiments of the present invention. As will be understood 
25 by those familiar with the art, the invention may be embodied in other 

specific forms without departing from the spirit or essential characteristics 
thereof. For example, the present invention may be implemented by a 
computer. Accordingly, the disclosure of the present invention is intended 
to be illustrative, but not Umiting, of the scope of the invention, which is set 
30 forth in the following claims. 



Claims 

I Claim: 

1. A method for approximating propagation delay through a logic 
device, the method comprising the steps of: 

5 (a) divide operation of the logic device into a first region and a second 

region, a boundary between the first region and the second being based on 
duration of input ramp to the logic device and amount of capacitive load 
driven by the logic device; 

(b) when the logic device operates in the first region, using a first 
10 formula to obtain a first value representing delay through the logic device, 

the first formula varying the first value based on the duration of the input 
ramp to the logic device and the capacitive load driven by the logic device- 
and; 

(c) when the logic device operates in the second region, using a second 
15 formula to obtain the first value, the second formula varying the first value 

based on the duration of the input ramp to the logic device and the capacitive 
load driven by the logic device. 

2. A method as in claim 1 wherein the first formula has a form as 
20 follows: 

DT = A + B*(INPUT RAMP) + D*(LOAD) 
where DT represents the first value, INPUT RAMP represents duration of 
the input ramp, LOAD represents capacitive load driven by the logic device, 
and A, B and D are coefficients determined using circuit simulations. 

25 

3. A method as in claim 2 wherein the second formula has a form as 
follows: 

DT = A+B*(CRITINRAMP)+D^(LOAD)+Z*(INPUT RAMP - CRITINRAMP) 
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wherein CRITINRAMP represents INPUT RAMP at the border between the 
first region and the second region for the LOAD, and Z is a determined 
using circuit simulations 

5 4. A method as in claim 1 wherein the first formula has a form as 

follows: 

DT = A + FAT*B*(INPUT RAMP) + D*(LOAD) 
wherein DT represents the first value, INPUT RAMP represents duration of 
the input ramp, LOAD represents capacitive load driven by the logic device, 
10 and A, B, D and PAT are coefficients determined using circuit simulations. 

5. A method as in claim 4 wherein the second formula has a form as 
follows: 

DT = DT(CRITINRAMP) + Z * (INPUT RAMP - CRITINRAMP) 
15 wherein CRITINRAMP represents INPUT RAMP at the border between the 
first region and the second region for the LOAD, Z is a coefficient 
, determined using circuit simulations, and DT(CRITINRAMP) is the first 
value when the logic ramp has an input ramp equal to CRITINRAMP and a 
load equal to LOAD* 

20 

6. A method as in claim 1 wherein the method includes calculating a 
next input ramp for a next logic device coupled to an output of the logic 
device, the method additionally comprising the steps of: 

(d) when the logic device operates in the first region, using a third 
25 formula to generate a second value representing duration from beginning of 
the input ramp of the logic device to an end of an output ramp of the logic 
device, the second value being based on the length of the input ramp to the 
logic device and the amount of the load driven by the logic device; and; 



(e) when the logic device operates in the second region, using a fourth 
formula to generate the second value; and, 

(f) subtracting the first value from the second value to obtain the next 
input ramp. 

5 

7. A method as in claim 6 wherein the first formula has a form as 
follows: 

DTo = Ao + B*(INPUT RAMP) + Dq*(LOAD) 
where DTq represents the first value, INPUT RAMP represents duration of 
10 the input ramp, LOAD represents capacitive load driven by the logic device, 
and Ao, B and Do are coefficients determined using circuit simulations. 

8. A method as in claim 7 wherein the third formula has a form as 
follows: 

15 DTi = Ai + B*(INPUT RAMP) + Di*(LOAD) 

where DTi represents the second value, and Ai and Di are coefficients 
determined using circuit simulations. 

9. A method as in claim 8 wherein the second formula has a form as 
20 follows: 

DTo = Ao+BKCRrriNRAMP)+Do*(^ - 

CRITINRAMP) 

wherein CRITINRAMP represents INPUT RAMP at the border between the 
first region and the second region for the LOAD, and Z is determined using 
25 circuit simulations. 

10. A method as in claim 9 wherein the fourth formula has a form as 
follows: 
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DTi * Ai + B*(INPUT RAMP) + Di*(LOAD). 

11. A method as in claim 6 wherein the first formula has a form as 
follows: 

5 DTo = Ao + FAT*B*(INPUT RAMP) + Do*(LOAD) 

wherein DTo represents the first value, INPUT RAMP represents duration 
of the input ramp, LOAD represents capacitive load driven by the logic 
device, and Ao, B, Do and FAT are coefficients determined using circuit 
simulations. 

10 

12. A method as in claim 11 wherein the third formula has a form as 
follows: 

DTi = Ai + FAT*B*(INPUT RAMP) + Di*(LOAD) 
wherein DTi represents the second value, and Ai and Di are coefficients 
15 determined using circuit simulations. 

13. A method as in claim 12 wherein the second formula has a form 
as follows: 

DTo - DTo(CRITINRAMP) + Z * (INPUT RAMP - CRITINRAMP) 
20 wherein CRITINRAMP represents INPUT RAMP at the border between the 
first region and the second region for the LOAD, Z is a coefficient 
determined using circuit simulations, and DT( CRITINRAMP) is the first 
value when the logic ramp has an input ramp equal to CRITINRAMP and a 
load equal to LOAD. 



25 



14. A method as in claim 13 wherein the fourth formula has a form 
as follows: 

DTi = DTi(CRITINRAMP) + B * (INPUT RAMP ^ CRITINRAMP). 



15. A method as in claim 6 wherein in step (a) the boundary between 
the first region and the second occurs where for each value of the capacitive 
load, the output ramp is one half complete when the input ramp is complete. 

16. A method as in claim 1 wherein in step (a) the boundary between 
the first region and the second occurs where for each value of the capacitive 
load, an output ramp for the logic device is one half complete when the input 
ramp is complete. 
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